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different topological sectors. 
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1. Introduction 



The chirally improved Dirac operator is a generalized Dirac operator which approximately 
obeys the Ginsparg-Wilson (GW) relation. Its construction and implementation for QCD has been 
discussed elsewhere M, 01. So far this operator was used for light hadron spectroscopy in quenched 
calculations by the Bern-Graz-Regensburg (BGR) collaboration [Q]. In those studies it was found 
that smearing the gauge links was important (using HYP smearing) as it resulted in better chiral 
properties for the operator, viz. the spectrum showed less deviation from the Ginsparg-Wilson 
circle compared to the unsmeared case. Therefore we decided to use one level of smearing in our 
studies, too. 

Usual HYP smearing is not well suited for use in HMC and therefore we implement the re- 
cently introduced stout smearing. It is tailor made for HMC and is differentiable. The main differ- 
ence between the stout smearing and other kinds of smearing is the way in which the general 3x3 
complex matrix is projected back to SU(3). Instead of the usual Gram-Schmidt orthogonalization 
one takes the traceless anti-hermitian part of the matrix (which is an element of the su(3) algebra) 
and then raises it to the group by the exponential map PJ]. 

We now simulate full QCD (with two light fermions) by implementing the HMC-algorithm 
for Dq\ and present here our first results. Details of the implementation and performance of the 
HMC-updating are discussed in another contribution to these proceedings 



2. Setting the scale 

Our Dirac operator follows ^ but uses stout smearing for the gauge configurations as part of 
the definition. We therefore had to reconstruct the the chirally improved Dirac operator appropriate 
for the combinations of gauge couplings and fermion masses used. The fermionic force in the 
HMC trajectory was calculated using the smeared chirally improved Dirac operator 

Since the stout-smearing is (to our knowledge) being used for the first time to construct the 
Dei, we find it instructive to look at the resulting difference between HYP and stout smearing. In 
Fig. [T] we plot different sized Wilson loops computed from HYP-smeared and stout- smeared links 
obtained from identical bare link gauge configurations. Actual values of the plaquette are quite 
sensitive to the smearing parameter for stout-smearing. We chose this parameter to be isotropic and 
one that maximizes the plaquette expectation value. Nevertheless we found that stout-smearing has 
less effect on the loop expectation values than HYP-smearing. Since smearing is a local process, 
the lattice spacing, derived from long distance behavior of correlators like loops, should not be 
affected by smearing. However, we expect less fluctuations from HYP smeared configurations. 

As a first indicator we set the scale on the lattice by using the Sommer parameter. (We may 
change to using hadron masses at a later stage.) We derive this from HYP smeared Wilson loops 
of various sizes. Computing the static potential by 



V(r) = -— "—log 
h - h 



W(r,T 2 ) 
W(r,Ti) 



(2.1) 



and assuming the usual form of the potential V(r) = or + ji+c/r, we obtain the parameters a, fl 
and c by fitting. Then the Sommer parameter is as usual given by r^Ffo) = 1.65, which in this 
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Figure 1: HYP smeared (l.h.s.) and stout smeared (r.h.s.) Wilson loops of different sizes. Note that for 
a (4 x 4) Wilson loop the expectation value of the HYP smeared loop is ~ 5 times larger than the stout 
smeared one. (Data from run on line 2 in Table 1.) 
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Figure 2: Sommer parameter (ro) using HYP smeared Wilson loops. 8 3 x 16 lattice with lattice spacing 
a = 0. 135 fm and bare quark mass m q = 0.05 in units of a. 

case translates to ro = ^(1.65 + c)/o. To set the scale explicitly in fermi we use ro = 0.5 fm. The 
sizes of the Wilson loops are such that ro lies inside the range of r used in the fits. 

To make sure that the lattice spacing obtained in this way (cf. Table |l|) does not depend on the 
smearing procedure, we repeated the calculation with stout smearing. That reproduced the results, 
albeit with slightly larger error bars. 

3. Results 

In Table [T] we summarize our run parameters. More details are given in [Qj. 
Next we look at the correlators of the two lightest mesons, the jt and the p mesons. Here we 
measure only point-to-point correlators defined by 

0,(0, = £ tr ( 75 D c l (x, t : 0, 0) y 5 D c / (0, : x, t)) , (3.1) 
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L 3 x T 


am 




res 




At 


steps 


acc. 


HMC 
time 


a[fm] 


8 3 x 16 


0.05 


5.4 


10 


-10 


0.015 


50 


~93% 


700 


0.114(3) 


8 3 x 16 


0.05 


5.3 


10 


10 


0.015 


50 


~91% 


700 


0.135(3) 


8 3 x 16 


0.08 


5.4 


10 


10 


0.015 


50 


~93% 


700 


0.138(3) 


12 3 x 24 


0.05 


5.3 


10" 


-8 


0.01 


100 


-87% 


200 


0.129(3) 



Table 1: Parameters, statistics and some results from our dynamical simulations. 



C P (0,t) = £ tr (n Del (x,t : 0, 0) y D c } {0, : x, t)) . (3.2) 

x,i= 1,2,3 

To extract the masses, the data were folded about the symmetry point and fitted to the func- 
tional form a cosh[m(? — T /2)]. For the 8 3 x 16 lattices, there seemed to be some contamination 
from the higher states and it was not clear if we were indeed observing the asymptotic decays for 
the correlators (cf. Fig. f|). On 12 3 x 24 lattices, such contaminations were less pronounced and 
within error bars we seem to obtain the asymptotic decays in the t-ranges shown in the Table ||. In 
this table we report on the masses for the 71 and the p meson as obtained from these correlators 
determined from fits in the given regions. The errors were obtained from a jackknife analysis. Due 
to the simple point-like sources the small distance behavior is dominated by excited states. This 
can be improved by using smeared sources. 

Since the Dq\ is only approximately a Ginsparg-Wilson operator, we have to a posteriori con- 
firm that the good chiral properties that we expect are indeed present in the operator. As a cross- 
check for that we compare the eigenvalue spectrum of Dqi for several configurations with dynami- 
cal fermions with those of a quenched simulation. The spectra have similar "fuzzyness" and are all 
close to the unit circle, as shown in [@]. 

For Dirac operators obeying the 75 hermiticity, real eigenmodes are the carriers of non- vanishing 
chirality and enumerate the total net topological charge. For massless exact GW fermions the 
lowest-lying eigenvalues are all exactly at the zero of the complex plane. For our approximate GW 
operator, even in the massless case, the lowest-lying eigenvalues are not exactly at zero but have a 
small real offset. Their total number defines the net topological charge. 

It is important to determine the amount of autocorrelation present in the system. This is ex- 
plored in |0|. Here we look at another closely related question, whether the algorithm samples all 
topological sectors of our model in the right way. We show a history of the number of close-to- 
zero modes (i.e., the absolute value of the topological charge) and symmetrized histogram of the 
topological charge (it is symmetrized simply because we did not compute the chirality of our real 
modes) in Fig. Q Data taken from different Markov chains are separated by dashed lines in the 
figure. We do observe tunneling (which seems to be a problem for the overlap operator) and a 
Gaussian-like shape of the distribution. 

4. Conclusions 

We have successfully implemented HMC with the chirally improved Dirac operator and our 
first results seem to indicate that the good chiral properties which were observed in the quenched 
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lattice-size a(fm) bare m q t - range am n y^jd.o.f am p % /d.o.f 



8 3 x 16 


0.135(3) 


0.05 


5-8 


0.520(15) 


~ 0.4 


0.808(14) 


~ 1.5 








6-8 


0.500(15) 


~0.05 


0.776(15) 


~0.1 


8 3 x 16 


0.114(3) 


0.05 


5-8 


0.598(58) 


~0.06 


0.881(26) 


~0.7 








6-8 


0.589(64) 


~0.02 


0.858(28) 


~0.06 


8 3 x 16 


0.138(3) 


0.08 


5-8 


0.621(18) 


~0.1 


0.860(13) 


~ 1.3 








6-8 


0.611(23) 


~0.01 


0.839(17) 


~0.1 


12 3 x 24 


0.129(3) 


0.05 


7-10 


0.384(21) 


~0.08 


0.575(37) 


~0.01 








8-10 


0.366(26) 


~ 0.04 


0.575(15) 


~0.06 



Table 2: n and p masses (in lattice units) from point to point correlators 

0.01 r , , , , , , , , a 




2 4 6 8 10 12 14 16 

t 

Figure 3: p (* , □) and % (A , x ) correlators on 8 3 x 16 lattices {a m q = 0.05). The upper line corresponds 
to lattice spacing a = 0.135 fm while the lower one is for a — 0.1 14 fm in both channels. 

case seem to be present in the dynamical case also. In contrast to the overlap operator, we do not 
have a strong barrier at the topological boundaries and we see quite often evidence of tunneling, 
manifested as "zero modes" of the Dqi- We also confirm (Fig. ||) that we have not undergone spatial 
deconfinement by looking at Wilson lines along any one of the shorter directions of the lattice. 

We estimate the masses of the 71 and p meson using point-to-point correlators. From quenched 
calculations ^ we know that 8 3 x 16 lattices (at a =0.15 fm) still have substantial size dependence 
and we expect the same to hold here too. However, the 12 3 x 24 lattices (~ 1.5 fm) should be better 
behaved. In this exploratory study we are at comparatively heavy quark masses m n /m p w 0.6. Runs 
using lighter masses are in progress. 
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Figure 4: History of real modes and symmetrized histogram of the topological charge. 
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Figure 5: Distribution of "spatial" Wilson lines for 8 3 x 16 lattice at a =0. 135 fm 
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